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Abstract. 

We describe the statistical properties of narrowband drifting auroral kilo- 
metric radiation ('striated' AKR) based on observations from the Cluster wide- 
band receiver during 2002-2005. We show that the observed characteristics, 
including frequency drift rate and direction, narrow bandwidth, observed in- 
tensity, and beaming angular sizes are all consistent with triggering by up- 
ward traveling ion solitary structures ('ion holes'). We calculate the expected 
perturbation of a horseshoe electron distribution function by an ion hole by 
integrating the resonance condition for a cyclotron maser instability (CMI) 
using the perturbed velocity distribution. We find that the CMI growth rate 
can be strongly enhanced as the horseshoe velocity distribution contracts in- 
side the passing ion hole, resulting in a power gain increase greater than 100 
dB. The gain curve is sharply peaked just above the R-mode cut-off frequency, 
with an effective bandwidth <50 Hz, consistent with the observed bandwidth 
of striated AKR emission. Ion holes are observed in situ in the acceleration 
region moving upward with spatial scales and speeds consistent with the ob- 
served bandwidth and slopes of SAKR bursts. Hence, we suggest that SAKR 
bursts are a remote sensor of ion holes and can be used to determine the fre- 
quency of occurrence, locations in the acceleration region, and lifetimes of 
these structures. 
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1. Introduction 

Auroral kilometric radiation (AKR) bursts exhibit a wide variety of fine structure as 
seen on frequency-time spectra. The cyclotron maser instability (CMI) [Wu and Lee, 
1979] is widely assumed to be the basic plasma mechanism responsible for the emission. 
This mechanism originally assumed a loss-cone electron velocity distribution function, but 
in situ observations in the acceleration region have shown that a 'horseshoe' or crescent 
distribution is more accurate [Louarn et al, 1990; Roux et al, 1993; DeLory et al, 1998; 
Ergun et al, 2000]. The horseshoe distribution, which arises naturally for electron beams 
in the presence of inhomogeneous magnetic fields [Speirs et al, 2005], provides a robust and 
efficient free energy source for the CMI mechanism, as shown both by model calculations 
[e.g. Pritchett, 1984; Pritchett et al, 1999; Bingham and Cairns, 2000] and in laboratory 
experiments [e.g., Speirs et al, 2005]. This mechanism has been applied not only to 
terrestrial AKR emission, but to many astrophysical environments in which the requisite 
conditions (low density, beamed electrons, inhomogeneous magnetic fields) are thought 
to be present, e.g. planetary magnetospheres [Zarka, 1998; Farrell et al, 2004], stellar 
magnetospheres [Louarn, Le Queau, and Roux, 1986; Kellett et al, 2002; Willes and Wu, 
2004], and even relativistic jets in active galaxies [Begelman et al, 2005]. 

Detailed physical models for the rich variety of observed spectral signatures are still 
an active area of research. AKR fine structure exhibits many morphologies on time- 
frequency spectra, including slowing drifting or nearly stationary features lasting tens of 
seconds [Gurnett et al, 1979], drifting features which may be interacting with each other 
[Pottelette et al, 2001], and periodically modulated or banded emission [Grabbe, 1982]. 
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Although some AKR bursts appear to be broadband emission extending over 100's of 
kilohertz [Hanasz et al, 2001], the total power is dominated by highly time- variable fine 
structures, especially during periods of enhanced bursts [Morozova et al, 2002], suggesting 
that perhaps all AKR radiation is a superposition of narrowband, short duration fine 
structures. 

Gurnett et al. [1979] first suggested that drifting AKR fine structure may be due to lo- 
calized sources rising and falling along auroral magnetic field lines, with emitted frequency 
equal to the local electron cyclotron frequency. Gurnett and Anderson [1981] suggested 
that the emission is triggered by electrostatic waves which drift along the field at the local 
ion-acoustic speed. An early detailed model which attempts to explain AKR fine structure 
is the tuned cavity model of Calvert [1982] in which the source region acts as a waveguide 
with sharp density boundaries. The source emits radiation in normal modes analogous to 
an optical laser, resulting in narrowband emission which drifts as the wave packet prop- 
agates into regions of varying width. A related idea is that of Farrell [1995] in which the 
cavity boundary is oscillating quasi-periodically. Under these conditions, broadband ra- 
dio photons will be stochastically accelerated, resulting in quasi-monochromatic discrete 
tones, whose frequencies slowly drift with changes in the cavity geometry. A difficulty 
with these models is that they require special conditions (e.g. sharp density boundaries, 
oscillating walls) which are not supported by observations [Pritchett et al, 2002]. 

Several authors have reproduced AKR fine structure using electromagnetic particle 
simulations to model the AKR source region. McKeen and Winglee [1991] used a one- 
dimensional electromagnetic particle-in-cell (PIC) code to model the cyclotron maser in- 
stability in a plasma with an inhomogeneous magnetic field. They found that radiation 
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is emitted in individual packets which combine to form drifting features that both rise 

and fall, consistent with some types of AKR fine structure. Both X- and O-mode drift- 
ing fine structures are created, consistent with observations [Benson et al, 1988] but not 
predicted in feedback models such as Calvert [1982]. A weakness of their model is that 
predicted drift rates are generally much higher than observed drifts. Pritchett et al. [2002] 
used a 2-dimensional PIC code with model parameters based on in situ FAST spacecraft 
measurements of the AKR source region. The simulated AKR is most strongly ampli- 
fied longitudinally and consists entirely of short-timescale fine structures. The predicted 
bandwidth (Au/lu ~ 10 -3 ) is somewhat larger than at least some observed AKR fine 
structure. 

More recently, Pottelette and colleagues have investigated the possible connection be- 
tween AKR fine structure and electron holes [Pottelette et al, 2001, 2003] and also tri-polar 
structures [Pottelette and Treumann, 2005]. AKR radiation from small packets ('elemen- 
tary radiation sources') associated with the holes are thought to be strongly amplified as 
the packets slow down and are reflected from the field-aligned electric field. Since the scale 
size of electron holes is very small (a few Debye lengths, Bounds et al. [1999]), the observed 
narrow bandwidth of AKR fine structures is easily accounted for. However, since neither 
the observed speeds or direction of electron holes (1000 - 2500 km s _1 always downward, 
[Bounds et al, 1999]) are representative of fine structure frequency drifts, there is not 
a straightforward relationship between the observed properties of electron holes and the 
dynamics of the AKR fine structure. 

In this paper, we present evidence that a particular type of AKR fine structure called 
striped or striated AKR [Menietti et al, 1996, 2000] is triggered by ion solitary structures 
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(ion holes). In section 2 we present new observations of SAKR bursts, including new 
bandwidth and angular beamwidth measurements. In section 3 we calculate the change in 
CMI gain for a density depleted region with a pre-existing horseshoe velocity distribution 
when perturbed by a passing ion hole. We compare the observed properties of SAKR 
with those calculated by radiation generated from the CMI instability in an ion hole. In 
section 4 we discuss some inferred properties of ion holes, such as lifetimes, probability 
of occurrence as a function of location, and radial dependence of speed, that cannot be 
inferred from single spacecraft in situ measurements. 

2. Observed properties of SAKR 

Striated AKR (SAKR) consists of trains of narrowband drifting bursts with negative 
slopes in the range -2 to -20 kHz s _1 , corresponding to upward-traveling sources with 
speeds between 100 - 1,000 kms -1 . Figure 1 shows a typical dynamical spectrum of SAKR 
bursts observed with the Cluster wideband (WBD) plasma wave instrument [Gurnett et 
al, 1997] on the four Cluster spacecraft. The bursts were observed in the 125-135 kHz 
band between 16:14:20 - 16:14:35 UT on 31 August 2002. Although the spectra are 
similar on all four spacecraft, there are clear differences in individual striations, indicative 
of angular beaming on the scale of the projected spacecraft separations, as shown in the 
inset at lower left. Also, individual bursts have very narrow bandwidths and are nearly 
always observed in groups with spacing between bursts of 30-300 ms. These characteristics 
are described in more detail below. 
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2.1. Occurrence frequency 

The occurrence probability of SAKR emission is quite low at frequencies above 100 
kHz: We detect SAKR bursts in less than 1% of all WBD spectra observed when the 
spacecraft was above 30°magnetic latitude (Note that below this latitude, there is often 
shadowing by the Earth's plasmasphere.) The occurrence probabilities were computed by 
dividing the number of dynamic spectra (length 52 sec) for which SAKR emission was 
clearly detected by the total number of spectra over a one year interval (July 2002 to 
2003). This overestimates the actual occurrence probability for SAKR since we do not 
correct for the fractional time within each spectrum that SAKR is present. 

There is a strong inverse correlation with observing frequency as shown Figure 2. There 
are only few detections in several hundred hours of data in the 500-510 kHz band (detection 
probability p < 0.01% ), while p ~ 0.2% in the 125-125 kHz band. Menietti et al. 
[1996, 2000] also studied SAKR emission using the plasma wave instrument (PWI) receiver 
on the Polar spacecraft. They also found an inverse correlation with observing frequency, 
with a 6% detection rate in the 90 kHz band (Fig. 2, shaded bar) and comparable 
probabilities and overlapping frequencies. By contrast, non-striated AKR emission is 
frequently detected, especially during sub-storm onset. The occurrence probability of 
AKR is nearly 30% for geomagnetic index K p in the range 1<K P < 3, with the highest 
occurrence frequency in the winter polar regions [Kumamoto et al, 1998]. 

2.2. Frequency drift and speed 

SAKR bursts have negatively sloped, nearly linear morphologies, at least within the 
normally sampled bandwidth (9.5 kHz) of the WBD receivers. The majority of the slopes 
range between -2 kHz s -1 and -8 kHz s _1 . Assuming the emission frequency is identified 
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with the local electron gyrofrequency, the speed of the stimulator as a function of frequency 
and slope can be written 



where V is the stimulator speed (kms^ 1 ), a is the observed slope of the SAKR burst 
(kHz s" 1 ), / is the observed frequency (kHz), and we have assumed the source moves 
upward along a dipolar magnetic field line at magnetic latitude \m = 70°. Figure 3 shows 
a histogram of the observed slopes of 650 SAKR burst events observed in the 125-135 
kHz band along with the derived trigger speeds (top x axis). The mean slope is -5.6 kHz 
s -1 , with more than 90% of the slopes in the range -2 to -8 kHz s _1 . The corresponding 
trigger speeds, using equation 1, are in the range 76 - 303 km s _1 , with a mean value 213 
km s" 1 . These results are similar to those reported by Menietti et al. [2000]. 

Figure 4 shows a group of SAKR bursts detected on one spacecraft while operating the 
WBD instrument in a wider bandwidth (77 kHz) mode. Over this wider bandwidth, the 
striations are not linear, but rather have a curved, frequency-dependent slope. The curve 
can be fit by assuming a source radiating at the local electron cyclotron frequency and 
moving with constant speed along a field line in a dipolar magnetic field. The overlaid 
white line shows the expected trace for a source moving upward at a constant velocity of 
300 kms -1 . SAKR bursts with similar frequency-dependent behavior are seen in Plate 1 
of [Menietti et al, 2000]. 

The altitude range of the SAKR locations are shown on the right y-axis. For this 
example, individual bursts originated near 6500 km altitude and moved upward at 300 
kms -1 to an altitude near 8100 km, implying a lifetime of several seconds. We have 
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examined several other SAKR events with wide bandwidth and have found that all exhibit 
similar lifetimes and have nearly constant speed, although the speed varies with epoch. 

2.3. Coeval broadband AKR emission 

SAKR bursts are almost always detected superposed on broader band AKR emission. 
This can be clearly seen in both Figure 1 and 4, as well as in Plate 1 of Menietti et al. 
[2000] . This appears to be a universal characteristic of SAKR, except in cases where the 
SAKR intensity is so low that underlying broadband emission may have been undetected. 
Background AKR emission coeval with SAKR bursts is consistent with the hypothesis 
that ion holes trigger the SAKR bursts (section 3.2) since the CMI gain in the source 
region is significant even in the absence of the ion holes. 

2.4. Bandwidth 

One of the most unusual features of SAKR bursts is their extraordinarily narrow band- 
width. The top panel of Figure 5 shows a single isolated SAKR burst. The middle panel 
shows the same feature, but after de-trending by subtraction of a 'chirp' signal of best- 
fit constant negative slope. The lower panel shows the full width at half maximum of 
a Gaussian fit along the frequency axis of the de-trended signal after summing in time 
intervals of 37 ms each. The resulting bandwidths, in the range 15 - 22 Hz, are much 
narrower than previously reported either observationally (e.g. Gurnett et al. 1979) or 
resulting from model calculation of AKR emission (e.g. Pritchett et al. 1999; Yoon and 
Weatherwax 1998) although Baumback and Calvert [1987] also reported AKR bandwidths 
as small as 5 Hz. (In the latter paper, the bursts do not appear to be SAKR emission). 
Other SAKR bursts we have examined have bandwidths ranging from 15 to 40 Hz. 
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The narrow bandwidth implies a small source extent along the z (B-field) direction, viz, 



Az = 0.55 km ■ 



100 kHz 



AL 

10Hz 



(2) 



/ 

where Az is the source extent along the B field and Af is the observed bandwidth. For 
/ = 130 kHz and Af = 20 Hz we obtain Az = 0.76 km. This is considerably smaller than 
the lateral extent of the AKR source region and indicates the trigger for SAKR must have 
a dimension along the magnetic field of order 1 km. 

2.5. Angular Beam Size 

The Cluster spacecraft array has the unique ability to simultaneously sample the flux 
density of individual AKR bursts at four widely separated points in space. By comparing 
the flux density on pairs of spacecraft, we can estimate the average angular beam size 
of individual bursts. The beaming probability is defined using the following algorithm. 
For each pair of spacecraft during a given observation, we calculate the projected angular 
separation between spacecraft as seen from a location situated above the magnetic pole 
of the hemisphere being observed, at a height corresponding to the electron cyclotron 
frequency at the center frequency of each observing band (e.g., 2.35 R e for the 125-135 
kHz band). This constitutes an average AKR location in that band without regard to 
location on the auroral oval. We next correct for differential propagation delay by shifting 
the waveform data from each spacecraft to the distance of the nearest spacecraft. We 
then divide the time-frequency spectrum from each spacecraft pair (52 sec duration for 
125-135 kHz band) into data 'cells' 19 msec x 53 Hz in size. This window was chosen to 
match the observed bandwidth of SAKR (cf. section 2.4) but the beaming results are not 



DRAFT 



February 5, 2008, 1:26pm 



DRAFT 



MUTEL ET AL.: STRIATED AKR: A REMOTE TRACER OF ION HOLES X - 11 

very sensitive to the data window size. We computed the angular beam size with data 
windows factors of two smaller and larger, and the results did not differ significantly. 

The intensities are normalized to correct for differing distances between individual space- 
craft and the AKR source. We then omit from further analysis all data cells whose in- 
tensities are below a threshold, arbitrarily chosen to be 10 dB below the maximum flux 
density for that spectrum. Finally, we compare intensities in each pair of data windows, 
assigning a weight 1 to pairs for which the intensities are within 10 dB of each other, and 
otherwise. The overall beaming probability for each angular separation interval is the 
sum of the beaming weights divided by all cell pairs. 

Note that this scheme is susceptible to overestimation of the beaming probability, since 
it is possible that independent AKR sources will illuminate separate spacecraft at the 
same frequency and time interval. This 'confusion' problem is smaller for SAKR emission 
since it is often clear from the burst morphology on a time-frequency spectrum that only 
one source contributes to a given data cell at one time, whereas with normal AKR there 
are very often several intersecting sources which contribute to a given data cell. 

In Figure 6 we plot the beaming probability versus angular separation for 651 SAKR 
bursts in the 125-135 kHz and 250-260 kHz bands. We have fitted a Gaussian function 
to the 125 kHz band observations using a least-squares fitting algorithm. The resulting 
full width at half maximum (FWHM) angular size is = 5.0° (solid angle Q = 0.006 sr). 
This is surprisingly small compared with most previously published observations of AKR 
beam size (e.g. Green and Gallagher [1985]), who reported beaming solid angles of 4.6 
sr and 3.3 sr at 178 and 100 KHz respectively. They made angular beaming estimates 
by comparing time-averaged spectra observed using two satellites (Hawkeye and IMP-6) 
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which simultaneously observed AKR bursts while the spacecraft were both over the same 
polar region. Since the time resolution used for the Hawkeye-IMP6 spectrum comparison 
(several minutes) far exceeds the time-scale of individual AKR bursts, their measured 
angular beam is actually a measure of the ensemble-averaged sky distribution of AKR 
bursts over a several minute time-scale rather than the angular beam size of individual 
AKR emission sources. This is the confusion problem mentioned above. 

An important unanswered question is the 2-dimensional structure of the SAKR burst 
angular beam pattern: It is asymmetric, or perhaps a hollow cone as suggested by Calvert 
[1987]? Since the spatial frequency coverage of the Cluster spacecraft array is often nearly 
one-dimensional at high magnetic latitude, it is difficult to analyze the 2-dimensional 
structure of individual bursts. The angular beaming probability plot shown in Figure 
6 includes measurements over a range of baseline orientations. However, a preliminary 
analysis of beaming probability grouped by baseline orientation did not reveal any obvious 
trends. We are presently analyzing a much a larger dataset consisting of a large variety of 
AKR emission and wider range of baseline orientations to investigate the 2-dimensional 
structure of the AKR angular beam. 

2.6. Individual SAKR Source Power Estimate 

The angular beaming observations, combined with measured flux density, allow a direct 
estimate of the average intrinsic power of individual SAKR bursts. For intense SAKR 
bursts, the observed square of the electric field intensity is 

I u ~ (1 - 10) x 1(T 12 (V/m) 2 Hz" 1 (3) 
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at a source-spacecraft distance d = 10R e . Converting to flux density, we obtain 

S„ = I u /Z ~ (0.3 - 3) x 10~ 14 W m" 2 Hz" 1 (4) 

where Z = 377 ohms is the impedance of free space. This flux density range is about one 
hundred times smaller than the flux density of intense AKR bursts reported by Benson 
and Fainberg [1991]. This is consistent with the conjecture that intense AKR bursts are 
the sum of many spatially distinct 'elementary radiation sources' [Pottelette et al, 2001] 
The power emitted at the AKR source is the isotropic power corrected by the angular 
beamsize of an individual SAKR burst 

P = And 2 -S v -Av— (5) 
An 

where Av ~ 20-50 Hz is the bandwidth of a burst, and Q ~ 0.006 sr, is the solid angle 
of the emission beam. Using these values, the resulting power is P ~ 1 - 10 W, much 
smaller than the previous estimates of P ~ 10 3 -10 4 W for a single 'elementary radiator' 
[Pottelette et al, 2001]. 

3. Connection with ion holes 

Ion solitary structures, also known as ion holes, are small-scale (~lkm) regions of 
negative electrostatic potential associated with upgoing ion beams. They are seen in 
spacecraft electric field measurements as symmetric bipolar parallel electric field structures 
with amplitudes 10 - 500 mV m -1 and timescales 3-10 ms. They were first detected in S3- 
3 spacecraft observations [Temerin et al, 1982] and have been subsequently been studied 
using in situ measurements in the acceleration region by several authors (e.g. Bounds et al. 
[1999]; Dombeck et al. [2001]; McFadden et al. [2003]). Ion holes travel upward at speeds 
between 75-300 kms -1 [Bounds et al, 1999; Dombeck et al, 2001] (although McFadden et 
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al. [2003] argue for somewhat higher speeds, in the range 550 - 1100 kms -1 ). The width 
of the waves increase with amplitude [Dombeck et al, 2001] which is inconsistent with 
small amplitude 1-d soliton models, but which supports a BGK-type generation mode 
[Muschietti et al, 2002]. 

3.1. Electron distribution function and CMI growth rate 

AKR radiation arises from wave growth resulting from the interaction of a radiation field 
with an electron velocity distribution having a positive slope in the direction perpendicular 
to the magnetic field (df/dvj_ > 0). The condition for waves of angular frequency uj and 
wave normal angle 6 to resonate with electrons with angular frequency uj cc and velocity v 
is given by [Wu and Lee, 1979] 



where uj ce is the electron cyclotron frequency, n is an integer, and 7 is the Lorentz factor 
of the relativistic electrons. The cyclotron maser instability (CMI) is the case n—\. If we 
assume that the electrons are mildly relativistic, as is observed in the AKR source region 
(E e ~ 10 keV, so 7 ~1.01), we can expand the Lorentz factor to obtain the equation for 
a circle in velocity space 



k\\vcos 6 — uj + n 








(6) 




(7) 



where the resonant circle's center is on the horizontal axis displaced by 




(8) 



and the radius of the resonant circle is 



2c 2 



UJ — UJ, 



ce 



2 



(9) 
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Analysis of recent FAST observations of AKR emission in the source region [Ergun et 
al, 1998, 2000; Pritchett et al, 1999] provides evidence, based on wave polarization, that 
the AKR k- vector direction at the source is nearly perpendicular to the B field (k\\ ~ 0). 
The E-field of the wave is polarized perpendicular to the ambient magnetic field which 
indicates that the wave is purely X-mode [Strange-way et al, 2001]. In the following, we 
explicitly assume k\\ = 0, so that the radius of the resonance circle becomes 



(10) 



where Sou = oo - oo ce . The growth rate of the CMI mechanism is given by calculating 
the imaginary part of the angular frequency, 



tUi = - — — f v\— — dv± (11) 



^ r vl di 

\.ujn e J dv± 

where the v± integral is performed on the closed circular path given by equation 4. 

The electron velocity distribution in the upward current acceleration region has been 
measured in situ by both Viking [Roux et al, 1993] and FAST [DeLory et al, 1998] 
spacecraft. It consists of an incomplete shell or 'horseshoe' shape in velocity space. The 
density of cold electrons (E << 1 keV) in this region is much smaller than the hot electron 
population which comprises the horseshoe component [Strangeway et al, 1998]. In this 
paper we assume that all electrons are in the hot component. 

We have modeled the observed velocity distribution using a simple analytic functional 
form 



f(v) = g{v) exp 



v - v N - 



(12) 



where A is the horseshoe width, f the horseshoe radius and the loss-cone function g(v) 

gM = l-/»-ed,( t ""^ ) (13) 
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where (3 is a dimensionless scaling factor and 9 is the characteristic opening angle of the 
loss cone. The model velocity distribution function is shown in Figure 7 along with a 
measured velocity distribution function from FAST [Ergun et al, 2000]. 

3.2. Triggering SAKR emission by an ion hole 

The range of possible interactions between background electrons and ion holes is very 
complex and varied [Eliasson and Shukla, 2005]. These include trapping of electrons be- 
tween ion holes, excitation and trapping of high-frequency Langmuir waves from electron 
streams resulting from ion hole collisions, and modification of the ion holes via the pon- 
deromotive force. In this paper we will tacitly assume that the ion holes are well-spaced 
and non-interacting, and that the only significant effect on the electron population is a 
transient speed decrease as the electrons traverse the hole's negative potential well. This 
assumption is motivated by the relatively simple structure of SAKR bursts and by the 
success in using this assumption to explain all significant observed properties of SAKR 
bursts. Other more complicated forms of AKR fine structure seen on dynamical spectra 
may well involve one or more of the complex interactions mentioned above. 

Since the electrons are magnetized (gyro-radius r ce ~50m, much smaller than the ion 
hole parallel scale size iiss ~1 km), the ion's hole's effect on the perpendicular velocity 
component is negligible (a small ExB drift, Yd ~ 40 m s _1 ). However, the parallel 
velocity component will experience a significant decrease as the electrons are repelled by 
the negative potential well of the hole. To determine the magnitude of the effect, we have 
modeled the observed electric field structure of an ion hole using the analytic form 




(14) 
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Figure 9a shows an observed E-field of an ion hole along with a plot of the model E-field 
with E = 500 mV m -1 , and r = 2 ms. 

We used this model to calculate the velocities of electrons as they traverse the hole for 
a variety of impact parameters and initial speeds. For initial parallel speeds of 20,000 
and 30,000 kms" 1 , the speed of an electron traversing an ion hole at a minimum distance 
of 1 km from the center is shown in Figure 9b. As expected, the parallel velocity briefly 
decreases, so that the horseshoe velocity distribution is 'squeezed' on the parallel (hori- 
zontal) axis, as shown in Figure 8b. The speed decrease is proportional to v~ l as expected 
by a simple energy conservation argument: For an electric field E and ion hole diameter 
Lj/i, the potential well of the hole is approximately $ ~E-Lj^. The electron loses kinetic 
energy as it traverses the potential well, so that energy conservation requires 

tfQmuf + =0 (15) 

Solving for the velocity change, 

* V „ = -^1 (16) 

m e v 

For an ion hole electric field E ~ 300 mV m" 1 , diameter Lqmi ~ 2km, and vq = 20,000 
kms" 1 , the expected speed decrease is 5v\\ ~ 5,000 kms" 1 , in good agreement with the 
exact calculation. 

In order to calculate the power gain in the perturbed region, we need an estimate of 
the convective growth length 

L c = ^ (17) 

where V g is the group velocity of the wave and uji is the growth rate. The group velocity is 
very sensitive to the ratio of plasma to gyro-frequency and the detailed velocity distribu- 
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tion function. As a first approximation, we have used the cold plasma dispersion relations 
to estimate the group velocity V g = du/dk for frequencies near the R-mode cut-off fre- 
quency. We find that Vg is of order 300 - 1000 km s _1 , similar to the estimate used by 
Omidi and Gurnett [1982]. 

For a maximum growth rate uJi ~ 5000 s _1 (at u ce = 2n x 125 kHz), the corresponding 
convective growth length is L c ~200m. Hence, there are approximately 5-10 e-folding 
lengths in the region of the ISS, assuming a perpendicular physical scale of order 1-2 
km. This results in a maximum power gain of e 10-20 ~ 43-87 dB. The power gain 
required to amplify the background radiation to observed levels of normal (wideband) 
AKR was estimated by Omidi and Gurnett [1982] to be e 20 , which is at the high end of 
our calculation. However, as discussed above, since SAKR emission is highly beamed, the 
required power is a factor of ~100 smaller, so the requisite power gain is closer to e 15 = 
65 dB. This is comfortably within the range predicted by an ion hole trigger. 

Figure 10 illustrates the growth rate and power gain both inside and exterior to an ion 
hole for a range of relativistic electron energies and plasma densities. Figure 10 (a, c) show 
the growth rate and power gain respectively for a horseshoe velocity distribution with a 
radius V/c = 0.15 (E = 5.7 keV) and n e = 0.1, 0.2 and 0.5 cm" 3 , while Figure 10(6, d) 
shows the growth rate and gain for n e = 0.2 cur 3 and a range of radii V/c = 0.10, 0.15, 
and 0.20 (E = 2.5, 5.7, 10.2 keV). These values are typical of the values of electron density 
and energy observed in the auroral density cavity [Strang ew ay et al, 1998]. There are 
several noteworthy features of these plots: 

1. The range of electron densities which result in substantial CMI gain is very limited. 
The maximum density is constrained by the condition that the R-mode cutoff frequency 
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is less than the electron cyclotron frequency. This can by expressed by the inequality 

fpeij^e) 
fee 

where f pe and f ce are the electron plasma and cyclotron frequencies respectively, and 7 is 
the Lorentz factor. This inequality can also be written 

n e< (J^Y (JL) cm -3 (19) 

V 9kHz / \ie ) 

where Eq is the rest-mass energy of the electron (511 keV). For example, for E = 5 keV 
and /= 125 kHz, we find n e < 0.9 cm~ 3 . On the other hand, for very low densities 
(n e << 0.1 cm -3 ), the CMI gain decreases dramatically (cf. Fig. 10a) 

2. The gain is a sharply peaked function of frequency. For example, for n e = 0.3 cm -3 , 
E = 5.7 keV (Fig. 10c), a dynamic range of 30 dB near the peak (approximately that 
observed on dynamic spectra of SAKR bursts) corresponds to a fractional bandwidth 
fif/f ~ 0.0005, or Sf ~ 50 Hz at / = 125 kHz, in good agreement with observed band- 
widths (section 2.4). These extremely narrow bandwidths are characteristic of CMI gain 
curves with partial-ring velocity distributions [Yoon and Weatherwax, 1998]. 

3. Even outside the ion hole there is substantial CMI gain under some favorable con- 
ditions (e.g. 80 dB for the middle plot of Figure lOd). This may be the explanation for 
the background AKR seen in Figure 4 and discussed in section 2.3. 

3.3. Conversion efficiency 

Given the rather small volume of an ion hole, it is reasonable to ask if there is sufficient 
free energy available from the resonant electrons to power the observed SAKR emission. 
As shown in section 2.5, the observed flux levels and angular beamwidths of SAKR emis- 
sion correspond to radiated powers P ra d ~ 1 — 10 W. Using a an ion hole scale size Lih ~ 1 
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km and speed ~ 300 km s -1 , and mean electron density and energy n e = 0.5 cm -3 , E e 
= 5 keV respectively, the ratio of radiated power to the electron kinetic energy traversing 
a hole per unit time is 

= (i - io)w ^ Q 0Q2 _ Q Q2 (20) 

Ptot n e E e Lih Vih 

This range of conversion efficiencies is similar to previous models of AKR emission 
[e.g., Pritchett et al, 1999] and indicates that while a single ion hole extracts relatively 
little energy from the ambient electron population, a train of hundreds of holes traversing 
the same region could account for a significant modification of the electron distribution 
function. 

4. Discussion 

If SAKR bursts are stimulated by upward traveling ion holes in the acceleration region 
of the magnetosphere, then observations of these bursts provides a new technique to study 
characteristics of ion holes such as lifetimes and relative number and speed as a function 
of altitude, that are impossible to measure with single in situ spacecraft . 

1. The 77 kHz bandwidth observations of SAKR bursts indicate that ion holes propa- 
gate upward for more than 1,000 km, implying lifetimes of a few seconds. This is much 
longer than estimates of solitary wave lifetimes derived from PIC simulations of solitary 
waves generated by the two-stream instability [Crumley et al, 2001], which have lifetimes 
r ~ (100 - 1500) Up' 1 -5-75 ms. 

2. The majority of ion hole speeds derived from SAKR observations are in the range 75 - 
400 km s" 1 (Fig. 3). These speeds are in very good agreement with in situ measurements 
of ion hole speeds from Polar (75 - 300 km s _1 ) at altitudes between 5500 km - 7000 
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km [Bounds et al, 1999; Dombeck et al, 2001]. However, McFadden et al. [2003] finds 
somewhat higher ion hole speeds (550-1100 km s^ 1 ) based on data from the FAST satellite 
at altitudes between 3000 km - 4000 km. It is possible that at speeds above V ~ 600 
km s _1 (slopes >20 kHz s _1 ) SAKR bursts, especially in closely spaced groups, would be 
undetected since they would merge into quasi-continuous emission on dynamic spectra. 

3. Fits to SAKR bursts assuming a constant speed source, (Fig. 4) indicate that 
ion holes propagate at nearly constant speed for their entire lifetime. Some numerical 
simulations [Crumley et al, 2001] show a significant change in ion hole speed as they 
evolve. This is not supported by our observations. 

4. SAKR bursts, and hence ion holes, are much more common at higher altitude, being 
more than one hundred times as common at 10,000 km (f ce ~ 90 kHz) than at 3,200 
km (f ce ~ 500 kHz) altitude, assuming ambient conditions favorable to the generation of 
SAKR are not dissimilar in this altitude range. 

5. Since SAKR bursts are almost always detected in groups with typical spacing 30 
- 300 ms., (section 2.6), this likely also applies to ion holes. Both laboratory plasma 
experiments [Frank et al, 2001] and observations of ion holes in the magnetosphere [e.g., 
Bounds et al, 1999] show trains of ion holes with spacings similar to the SAKR bursts. 

6. Finally, the uniformity of SAKR intensity and bandwidth over a large frequency 
range (Fig.4) implies that there is little evolution of the electric field intensity or spatial 
structure of ion holes over their lifetime. 
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5. Summary and Conclusions 

This paper summarizes the observed properties of a distinct form of auroral kilometric 
radiation fine structure called striated AKR, first described by Menietti et al. [1996]. We 
present new observational results using the WBD instrument on Cluster which characterize 
the bandwidth and angular beamsize of individual SAKR sources, as well as derived 
properties such as intrinsic power and speed along the magnetic field. 

Assuming the frequency of SAKR bursts can be identified with the local electron cy- 
clotron frequency, the speed and direction of SAKR sources calculated from their observed 
drift rates are very similar to those observed for ion solitary structures (ion holes) in the 
upward current region of the magnetosphere. Hence, we investigated whether SAKR 
bursts could be the result of enhancement of the cyclotron maser instability by ion holes. 
Using observed electric field signatures of ion holes in this region, we calculate the pertur- 
bation caused by the passage of an ion hole on a 'horseshoe' electron velocity distribution 
in a dilute (n e << 1 cm -3 ) plasma. The cyclotron maser instability is strongly enhanced 
inside the ion hole, with power gain exceeding 100 dB in a narrow frequency range just 
above the x-mode cutoff frequency. These characteristics are in excellent agreement with 
the observed bandwidth, speed, direction, and flux density of SAKR bursts. Alternative 
suggestions involving other types of solitary structures to explain AKR fine structure, such 
as electron holes [Pottelette et al, 2001] or tri-polar structures [Pottelette and Treumann, 
2005], are not consistent with SAKR properties, although they may be important in other 
types of AKR fine structures. 

If SAKR bursts are in fact triggered by ion holes, a number of derived properties of 
ion holes can be deduced which would be difficult or impossible to obtain using in situ 
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measurements. These include average lifetimes (a few seconds), evolution of propagation 
speed (nearly constant over lifetime of hole), and relative numbers versus location (much 
more common high in the acceleration region than near the base). 
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Figure 1. Frequency-time spectra of striated AKR bursts observed on 4 Cluster 
spacecraft on August 31, 2002 from 16:14:20 - 16:14:35 UT. The inset shows the projected 
angular separations of the spacecraft as seen from the source. 
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Histogram of AKR Rain Observations 
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Figure 2. Frequency of occurrence of SAKR as a function of observing frequency. The 



data for 90 kHz (hatched bar) is from Menietti et al. [2000]. 
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Figure 3. Histogram of observed SAKR burst slopes in the 125-135 kHz band (bottom 
x-axis) and derived trigger speed (top x-axis, from equation 1). 
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Figure 5. SAKR bandwidth, (a) SAKR burst observed with Cluster WBD, (b) Same 
burst, but detrended, (c) Measured bandwidth in 37 ms intervals. 
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Figure 6. Angular beamsize for 651 individual SAKR burst sources measured at 125 
kHz (red dots) and 250 kHz (green dots) using observed correlations between multiple 
spacecraft. The size of each point is proportional the number of bursts for that point. 
The fitted FWHM beamsize is 5.0°. 
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Figure 7. (a) Electron distribution function measured by FAST in AKR source region 



[Ergun et ai, 2000]. (b) Model distribution. 



DRAFT 



February 5, 2008, 1:26pm 



DRAFT 



X - 36 MUTEL ET AL.: STRIATED AKR: A REMOTE TRACER OF ION HOLES 




Figure 8. Top panel: horseshoe electron velocity distribution function outside an 
ion hole (a), and inside an ion hole (b). Black circle is the CMI resonant circle. The 
v || direction is the horizontal axis, Earthward to left. Bottom panel: weighted partial 
derivative df/dvj_ outside an ion hole (c), and inside an ion hole (d). The arrows indicate 
the parts of the growth line integral (equation 9) enhanced inside the hole. 
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Figure 9. (a) Observed electric field vs. time (solid line) as an ion hole passes by FAST 
spacecraft [McFadden et al, 2003], and best-fit model (dashed line) using equation (14). 
The inset at lower right shows a simple spherical model of an ion hole with inward electric 
field, (b) Calculated parallel velocity vs. distance from an ion hole for an electron passing 
through an ion hole with an impact parameter 1 km and initial velocities 20,000 km s _1 
(dotted line) and 30,000 kms" 1 (solid line). The assumed hole speed is 300 kms -1 . 
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Figure 10. (a) Temporal growth rate as a function of fractional frequency for electron 
horseshoe distribution radius V c /c = 0.15 (E = 5.7 keV) and electron density n e = 0.1 
cm _3 (A), 0.2 cm -3 (□), and 0.5 cm~ 3 (o). Solid lines are inside ion hole; dashed lines 
are outside ion hole. (b). Same as (a), but with fixed n e = 0.2 cm -3 and V c /c = 0.10 
(2.5 keV,D), 0.15 (5.7 keV,o), 0.20 (10.2 keV, A), (c) CMI power gain vs. fractional 
frequency, same parameters as (a), (d) CMI power gain vs. fractional frequency, same 
parameters as (c). 
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